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Abstract 

In general, the thermal-convection (Navier-Stokes and 
energy) equations are nonlinear partial differential equations 
that in most cases require the use of complex algorithms in 
combination with opportune discretization techniques for 
obtaining reliable numerical solutions. There are some cases, 
however, in which such equations admit analytical solutions. 
Such exact solutions have enjoyed a widespread use in the 
literature as basic states for determining the linear stability 
limits in some idealized situations. This review article 
provides a synthetic review of such analytical expressions 
for a variety of situations of interest in materials science 
(especially crystal growth and related disciplines), including 
thermogravitational (buoyancy), thermocapillary 
(Marangoni), thermovibrational convection as well as 
"mixed" cases and flow controlled via static and uniform 
magnetic fields. 
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Introduction 

Finding solutions to the Navier-Stokes equations is 
extremely challenging. In fact, only a handful of exact 
solutions are known. For the cases of interest in the 
present article (buoyancy, Marangoni, vibrational and 
magnetic convection), in particular, these exact 
solutions exist when the considered system is 
infinitely extended along the direction of the imposed 
temperature gradient or the orthogonal direction; in 
general, such solutions are regarded as reasonable 
approximations in the steady state of the flow 
occurring in the core of real configurations which are 
sufficiently elongated (the core is the region 
sufficiently away from the end regions, where the 
fluid turns around, to be considered not influenced by 
such edge effects). Even though limited to cases of 
great simpicity, these solutions have proven able to 
yield directly or indirectly insights and understanding 
that would have been difficult to obtain otherwise. 

Hereafter, first the attention is focused on the case of 
systems subject to horizontal heating (or to 



temperature gradient along the direction in which 
they are assumed to be infinitely extended), then other 
possible variants and configurations are considered in 
the section on "inclined systems" and subsequent ones. 

Analytic Solutions for Thermogravitational 
and Marangoni Flows 
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FIG. 1 SKETCH OF LAYER OF INFINITE EXTENT SUBJECT TO 
HORIZONTAL HEATING 

In the following three subsections (focusing on the 
Hadley flow, Marangoni flow and related hybrid 
states) the horizontal boundaries are assumed to be 
located at y=-l/2 and 1/2, respectively. The velocity 
components along y and z are zero (v=w=0) and the 
component along x is solely a function of y, i.e. u=u(y) 
(which means the assumption of "plane-parallel flow" 
is considered). Temperature depends on x and y (i.e. 
these solutions are essentially twodimensional). 
Rayleigh and Marangoni numbers are defined as 
Ra=GrPr=gPTyd 4 /va and Ma= RePr=aTyd 2 /|ia, 
respectively (where y is a rate of uniform temperature 
increase along the x axis, d is the distance between the 
boundaries, a is the thermal diffusivity, v is the 
kinematic viscosity, Pt is the thermal expansion 
coefficient, Pr the ratio between v and a, Gr and Re the 
Grashof and Reynolds numbers, respectively). 
Referring velocity and temperature to the scales a/d 
and yd, respectively and scaling all distances on d, the 
governing equations in nondimensional form 
(incompressible form) read 



V-V = 



cT 



= = -yp-y-[yy]+Prv 2 y-PrRtf r i g 
+ y_-[vr] = v 2 r 



(1) 
(2) 
(3) 
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where i g is the unit vector along the direction of 
gravity and the Boussinesq approximation has been 
used for the buoyancy production term in the 
momentum equation. 

In the presence of a free liquid/gas interface such 
equations must be considered together with the 
socalled Marangoni boundary condition, which 
neglecting viscous stress in the gas (here the dynamic 
viscosity of the gas surrounding the free liquid surface 
will be assumed to be negligible with respect to the 
viscosity of the considered liquid) reads 



3i 



-Ma V S T 



(4) 



where n is the direction locally perpendicular to the 
considered elementary portion of the free interface, v s 
the derivative tangential to the interface and V s the 
surface velocity vector. 

For the case considered here (steady flow with u=u(y) 
and v=w=0, as shown in Fig. 1) eqs. (l)-(3) reduce to : 



dp „ 5 « 
— = Pr — - 

dx dy 2 



dp 



= PtRaT 



dy 

dT d 2 T d 2 T 

u — = — - + — - 
dx dx dy 



(5) 
(6) 
(7) 



while eq. (4) can be rewritten as: 



(8) 



du , , dT 

— = -Ma — 
dy dx 

Assuming solutions in the general form (with Ma=0 or 
Ra=0 for pure buoyancy or Marangoni flows, 
respectively): 





u 




y = 


V 






w 





Ra gl (y) + Ma g 2 (y) 





T = x + Ra f,{y) + Ma f 2 (y) 



(9) 



(10) 



fx, fi, gi and g2 can be determined as polynomial 
expressions with constant coefficients satisfying 
specific equations obtained by substituting eqs. (9) and 
(10) into eq. (7) and into the equation resulting from 

cross-differentiation of eqs. (5) and (6) (i.e. <*^_ - dT ' 

dy 3 



Ra—)- 

dx 



Such system of equations read: 

d 2 L 



8 



dy 2 



(11a) 



Ra\ 



d 2 f 2 
' dy 2 
d 3 8l 

dy 3 



-Ma 



d 3 g 2 
dy 3 



= 



(lib) 
(12) 



to be supplemented with the proper boundary 
conditions: 



Kinematic conditions: 

solid boundary: u=0 — > gi(y)=g2(y)=0 (13) 
free surface: eq. (8) — > dgi 



^ and*i 

dy dy 



-1 



(14) 



Thermal conditions: 



Adiabatic boundary: ^L = ^^± = ^± = (15) 

dy dy dy 

Conducting boundary: T=x— > /i=/ 2 =0 (16) 

In addition, it must be taken into account that since the 
considered parallel flow is intended to model a slot 
with distant end walls, continuity requires that the net 
flux of fluid at any cross section of the slot be zero), i.e.: 



/ z / z / z 

judy = 0^> jg t dy= jg 2 dy = 



(17) 



Before going further with the description of solutions 
satisfying the system of equations (11-17), it is worth 
noting that some insights into their expected 
polynomial order can be given immediately on the 
basis of eq. (12). According to such equation, in fact, 
d 3 gi/dy 3 =l for Ma=0, while d 3 g2/dy 3 =0 for Ra=0, which 
leads to the general conclusion that the polynominal 
expression for u will be of the third order for pure 
buoyancy flow and of the second order for pure 
Marangoni flow, while (on the basis of eqs. (11)) the 
respective temperature profiles are of the fifth and 
forth orders in y. 

Thermogravitational Convection: The Hadley flow 

In line with the considerations above, in the case of 
liquid confined between two horizontal infinite walls 
with perpendicular gravity, eqs. (5)-(7) admit as an 
exact solution the following velocity profile: 



Ra 



u = ■ 



(18) 



generally known as the Hadley flow (Hadley, 1735) as 
it was originally used as a model of atmospheric 
circulation (see, e.g., Lappa, 2012) between the poles 
and equator (from a technological point of view, this 
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solution is also strongly relevant to the manufacture of 
bulk semiconductor crystals). 

The corresponding temperature profile changes 
according to the type of boundary conditions 
considered. For adiabatic walls it reads: 



Maf 2 1 
■Xl 3y +y ~4 



(20) 



T , RCl { 4 5 2 5 

T = x-\ v y y H 

120 { 6 16 



whereas for conducting boundaries it becomes: 



T = x H y y y H 

120 I 6 48 



(19a) 



(19b) 



The polynomial expressions / and g for such solutions 
are plotted in Figs. 2. 
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FIGS. 2 EXACT SOLUTION FOR THE CASE OF 
THERMOGRAVITATIONAL CONVECTION IN INFINITE LAYER 
WITH TOP AND BOTTOM SOLID WALLS: (a) VELOCITY 
PROFILE g(y); (b) TEMPERATURE PROFILE /(y) FOR ADIABATIC 
(SOLID) AND CONDUCTING (DASHED) BOUNDARIES. 

Marangoni Flow 

In the absence of gravity and replacing the upper solid 
wall with a liquid-gas interface supporting the 
development of surface-tension-driven (Marangoni) 
convection, the exact solution reads (Birikh, 1966): 



For adiabatic interface and insulated bottom wall the 
temperature distribution can be expressed as: 



Maf 
= x 



3>' 4 +2y 3 -V4;y- 



16 



(21a) 



48 V" " 2" T 
that for conducting boundaries must be replaced with: 

M "f 3v 4 +2v 3_3 2 _\_ 3^ (21b) 
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FIGS. 3 EXACT SOLUTION FOR THE CASE OF 
THERMOCAPILLARY CONVECTION IN INFINITE LAYER WITH 
BOTTOM SOLID WALL AND UPPER FREE SURFACE: (a) 
VELOCITY PROFILE g(y); (b) TEMPERATURE PROFILE f(y) FOR 
ADIABATIC (SOLID) AND CONDUCTING (DASHED) 
BOUNDARIES. 

The polynomial expressions / and g for such solutions 
are plotted in Figs. 3. 

Hybrid Buoyancy-Marangoni States 

As shown by eqs. (5)-(7), under the considered 
conditions (u^O, v=w=0), the nonlinear convective term 
of the momentum equation becomes zero and the 
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energy equation reduces to u=9 2 T/9y 2 , i.e. the 
governing equations are linear; as a consequence, 
more complex solutions can be built as superposition 
(addition) of other simple existing solutions. 
Along these lines, in the presence of vertical gravity 
and a liquid-gas interface supporting the development 
of surface-tension-driven (Marangoni) convection, the 
velocity profile can be obtained as the sum of two 
terms, the first term corresponding to the pure 
buoyancy-driven flow and the second term 
representing the contribution of a pure 
thermocapillary-driven flow, which provides a simple 
theoretical explanation for the general form given by 
eqs. (9) and (10). 

The resulting shear flow is set up by a combined effect 
of buoyancy and viscous surface stress due to the 
temperature dependence of surface tension. The first 
contribution related to pure buoyancy flow, however, 
exhibits some differences with respect to eq. (18) as the 
upper solid wall considered earlier must be replaced 
with a stress-free boundary (see Fig. 4); the velocity 
profile reads: 

u=—(sy 3 -3y 2 -3y+ l 
48 I 



(22) 



Hence, the resulting profile for mixed gravitational- 
Marangoni convection has the form: 



Ra 
48 



8y 3 -3y -3y + 



1 



Ma ( 



3y 2 + y 



(23) 



The associated temperature distribution is: 
1 (Ra 



T = x + - 



48 20 



8 /_ 5 /_ 10 / + V +5 W_ 3 ( 1 
I ib \ 2 



-Ma 



3y 4 + 2y 3 -V- 



3 5(1 

2 ), -i6 + r + 2 



(24) 



The cases of adiabatic and conducting horizontal 
surfaces can be obtained by setting in this relation r)=0 
and r)=l, respectively. 

Equation (24) can be extended to the more general case 
in which the bottom is conducting and a Biot number 
is introduced to describe heat transfer at the top free 
surface by assuming r|=Bi/(l+Bi). The solutions for 
symmetrical (top and bottom) conditions can again be 
recovered as limiting cases when Bi— »0 and Bi— ho: 
When Bi— >oo, in fact, r|— >1 and eq. (24) reduces to the 
conducting case; when Bi-»0 (r|-»0), it reduces to the 
insulating boundary conditions for the temperature 
profile. 
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FIGS. 4 EXACT SOLUTION FOR THE CASE OF PURE 
THERMOGRAVITATIONAL CONVECTION IN INFINITE LAYER 
WITH BOTTOM SOLID WALL AND UPPER STRESS-FREE 
SURFACE: a) VELOCITY PROFILE g(y); b) TEMPERATURE 
PROFILE /(y) FOR ADIABATIC (SOLID) AND CONDUCTING 
(DASHED) BOUNDARIES. 

General Properties 

The present section is devoted to illustrate some 
fundamental properties of the solutions introduced in 
the preceding sections with respect to the well-known 
Rayleigh's condition, which so much attention has 
enjoyed in the literature as a means for gaining 
insights into the fluid-dynamic instabilities of parallel 
flows in the limit as Pr^O (Rayleigh, 1880; Lin, 1944; 
Rosenbluth and Simon, 1964; Drazin and Howard, 
1966). 

Let us recall that such a theorem reads: In a shear flow a 
necessary condition for instability is that there must be a 
point of inflection in the velocity profile u=u(y), i.e. a point 
where d z u/dy 2 =0. 

The second derivative of solution (18) gives: 



4 



Review of Applied Physics (RAP) Volume 1 Issue 1, December 2012 



www.seipub.org/rap 



Ra y 



(25) 



d 2 u 
df 

that means the profile has an inflection point at 
midheight (y=0) and satisfies the Rayleigh's necessary 
condition. 

For pure buoyancy flow with upper free surface, i.e. 
solution (22), the second derivative reads: 
Ra 



d 2 u 



dy 2 48 



(48 y- 6) 



(26) 



that gives the inflection point at a sligtly different 
position (y=l/8). 

The inflection point, however, is no longer present for 
the case of pure Marangoni flow since: 
d u 3 



dy 2 



-Ma * 



(27) 



that means Marangoni flow solutions of the type given 
by eq. (20) do not satisfy the Rayleigh's necessary 
condition. 

The most interesting case in this regard is, perhaps, 
given by the mixed state represented by eq. (23). In 
such a case: 
d 2 u Ra 
^"48 



(48j-6-72W) 



(28) 



that makes the location of the inflection point a linear 
function of the nondimensional parameter W=Ma/Ra: 
12W+1 

y = — o — 



(29) 
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FIG. 5 VELOCITY PROFILE gi(y)+Wg2(y) FOR THE CASE OF 
MIXED THERMOGRAVITATIONAL-THERMOCAPILLARY 
CONVECTION IN INFINITE LAYER WITH BOTTOM SOLID 
WALL AND UPPER STRESS-FREE SURFACE: THE INFLECTION 
POINT DISAPPEARS FOR W=0.25. 

The inflection point disappears when y>l/2, i.e. for 
W>l/4. Accordingly, these solutions satisfy the 



Rayleigh's necessary condition only if 0<W<l/4 (see 
Fig. 5). 

The Infinitely Long Liquid Bridge 

For pure Marangoni flow, the equations of thermal- 
convection admit exact solutions also for the case of a 
liquid bridge assumed to be infinitely extended along 
the axial direction z (Xu and Davis, 1983). For such a 
case the solution is axisymmetric and reads: 



w - 



Ma 
2 



r 



T = -z- 



Ma 
32 



(30) 



(31) 



The interface (adiabatic) is assumed to be located at 
r=l. The Marangoni number is defined as 
Ma=RePr=aTyR 2 /|ia (where y is the rate of constant 
temperature increase along the z axis, R is the radius 
of the liquid bridge). Velocity and temperature are 
referred to the scales a/R and yR, respectively 
(moreover, all distances are scaled on R). 

According to eq. (30), the velocity component along 
the axis of the liquid bridge is solely a function of z 
(w=w(z)). Temperature is a function of both axial and 
radial coordinates (z and r). 

Inclined Systems 

The solutions given before for the case of a horizontal 
layer of fluid with no-slip walls subjected to a uniform 
temperature gradient along x can be extended to the 
more general case in which the layer is inclined by an 
angle 90°-9 with respect to the horizontal direction 
(see Fig. 6). 

y 
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Hot 



FIG. 6 SKETCH OF A LAYER OF INFINITE EXTENT INCLINED 
WITH RESPECT TO THE HORIZONTAL DIRECTION 
(TEMPERATURE GRADIENT ALONG THE x DIRECTION; 
ANGLE BETWEEN VT AND g). 

In particular, the analytical solution for such a case has 
distinct expressions depending on the projection of the 
temperature gradient to the gravity vector, i.e. 
according to whether from a global point of view the 
layer tends to behave as a system heated from below 
or from above (depending on the sign of 90°-8 in Fig. 
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6), and depending on the type of thermal conditions 
along the walls. 

In the following these expressions are given first for 
the case of adiabatic walls (eqs. (32)-(35)); then the 
configuration with conducting boundaries is 
considered (eqs. (38)-(41)). 



Heating from below (0°<9<90°) 



Rasm(0) 



16 



sinh(2£y) sin(£) - sinh(£) sin(2£y) 



' 4 2 (sinh(^) cos(^) + cosh(^) sin(£))_ 



(32) 



jt + — tan(0) 
2 



2y 



sinh(2gy) sin(^) + sinh(^) sin(2gy) 
sinh(^) cos(^) + cosh(^) sin(£) 



(33) 



Heating from above (90°<9<180°) 



u(y) ■■ 



Ra sin(#) 
16 



cosh(24y) sin(2£y) sinh(^) cos(^) 
? (sinh(£) cosh(£) + sin(£) cos(£)) " 



sinh(2£y) cos(2£y) cosh(^) sin(£) 
? (sinh(|) cosh(£) + sin(£) cos(£))_ 



(34) 



x + -tan(&) 



2y 



sinh(2£y) cos(2^y) sinh(^) cos(^) 



^(sinh(^) cosh(^) + sin(£) cos(£)) 
cosh(2gy) sin(2gy) cosh(^) sin(|) | (35) 
4{smh(4) cosh(^) + sin(£) cos(#)) 

where 



£ = |[/?acos(0)] I/4 for9<90° 

$ = ^rA- Racos(0)] U4 for 9>90° 
2V2 



(36a) 
(36b) 



Interestingly, some fundamental insights into this kind 
of solutions can be obtained (as noticed by Delgado- 
Buscalioni and Crespo del Arco, 2001) considering the 
vorticity-production term of the vorticity balance 
equation. For the present case such term reads: 



Ra 



dT . dT 

— sin(<9) + — cos(r?) 

dx dy 



(37) 



The evolution of these flow profiles with Ra is ruled 
by the balance of dissipation and production of 
vorticity by buoyant forces. According to eq. (37) 
above, the production of vorticity due to the y and x 
components of buoyancy are proportional to Ra 
sin(9)3T/3x and Ra cos(9)3T/5y, respectively. 

The temperature y-gradient is created by the flow 
advection and at low enough values of Ra it is 
negligibly small; therefore, at small Ra and for any 
(not vertical) inclination the flow is generated solely 



by the y component of buoyancy (at low values of Ra, 
in the conducting regime the cross-stream temperature 
gradient is vanishingly small and the vorticity is 
generated by the cross-stream component of gravity, 
at a rate given by sin(9)3T/3x; this induces a cellular 
flow whose y-dependence coincides for 9=90° with the 
profile (6)). 

As Ra increases, the streamwise advection creates an 
increasing (positive) temperature gradient along the y 
axis, which acts as another source of motion owing to 
the presence of the streamwise component of 
buoyancy; as explained before, this term produces 
vorticity at a rate given by cos(9)3T/3y and hence its 
effect greatly depends on the range of the inclination 
angle. 

When heating from above (9>90°) the effect of the axial 
buoyancy is to suppress the convection in the center 
part of the layer, as long as cos(9)3T/5y>0 while 
sin(9)5T/3x<0. For large enough Ra and 9>90°, the flow 
is confined to small regions near the walls where 
3T/5ysO. 

On the contrary, if the cavity is heated from below 
(9<90°), both sources of vorticity have the same sign 
and as Ra increases a positive feedback loop between 
u(y) and T(y) occurs: any increment of the flow 
intensity increases the cross-stream temperature 
gradient, which in turn, enhances the intensity of the 
flow. 

For conducting boundaries the expressions for heating 
from below and from above read: 

Heating from below (0°<9<90°) 



u = 



Rasin(0) 



16 



sinh(2fo) sin(£) - sinh(^) sin(2fo) 1 
£ 2 (sinh(£)sin(<r)) ' ' \ 



T = x + -tm(0) 



2y 



(38) 

1 sinh(2£y)sin(£) + sinh(£)sin(2&0 



sinh(^) sin(£) 



(39) 



Heating from above (90°<9<180°) 



u(y) = 



Ra sin(0) 



16 



(40) 



cosh(2£y) sin(2£y) sinh(^) cos(^) 
4 2 (sinh 2 (4) cos 2 (4) + cosh 2 (£) sin 2 (£)) 

sinh(2£y) cos(2£y) cosh(^) sin(£) 
4 2 (sinh 2 (4) cos 2 (4) + cosh 2 (4) sin 2 (4))_ 

sinh(2£y) cos(2£y) sinh(^) cos(^) 
y (sinh 2 (4) cos 2 (4) + cosh 2 (4) sin 2 (£)) 

cosh(2£y) sin(2^y) cosh(^) sin(£) 
(sinh 2 (4) cos 2 (4) + cosh 2 (4) sin 2 (4))_ 

with t, given by eqs. (36). 



T = x+-tan(<9) 



(41) 
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FIG. 7 SKETCH OF A LAYER OF INFINITE EXTENT INCLINED 
WITH RESPECT TO THE HORIZONTAL DIRECTION WITH 
HEATING APPLIED THROUGH THE BOTTOM WALL 
(TEMPERATURE GRADIENT IMPOSED ALONG THE y 
DIRECTION). 

Equations (32)-(35) and (38)-(41) provide analytical 
solutions if the imposed AT is parallel to the walls as 
shown in Fig. 6; thermo gravitational convection in an 
inclined layer, however, also admits exact solutions if 
the imposed T is primarily perpendicular to the 
walls i.e. if such temperature gradient acts across the 
thickness of the layer (i.e. layer with upper and lower 
walls kept at uniform different temperatures as shown 
in Fig. 7). Such a solution reads: 



Ra 



sin(6>) 



y 



y 

4y 



(42) 



with Ra=gPATd 3 /va. 



The corresponding distribution of temperature along y 
is governed by diffusion only and is approximately 
uniform throughout the plane of the layer. 

Velocity (u) tends to zero (as expected) in the limit as 
6— »0 (in such a case, in fact, it is a well-known fact that 
convection arises only if a given threshold of the 
Rayleigh number is exceeded). 
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FIG. 8 SKETCH OF A TRANSVERSELY HEATED VERTICAL 
CAVITY IN THE LIMIT AS THE VERTICAL SCALE OF MOTION 
TENDS TO INFINITY. 

It is also worth noticing that for 6=90° one recovers the 
idealized case of a transversely heated vertical cavity 
(in the limit as the vertical scale of motion tends to 
infinity, see Fig. 8): 



v = 



V 



4. 



(43) 



this profile is generally referred to as "conduction- 
regime" solution (owing to the associated temperature 
profile along x, that as mentioned before is linear). 

Convection with Vibrations 

Pure Thermovibrational Flow 

Let us recall that thermovibrational convection can be 
regarded as a "variant" of the standard 
thermogravitational convection for which the steady 
Earth gravity acceleration is replaced by an 
acceleration oscillating in time with a given frequency 
(Lappa, 2004, 2010). 

Disturbances induced in a fluid by a sinusoidal 
displacement of the related container along a given 
direction ( „ is the related unit vector) 



s(t) = bsin(&> t)n 



(44) 



where b is the amplitude and ©=2nf (f is the frequency) 
induce an acceleration: 



g{t)=g„ sin (« 



(45) 



where g = b afn ', which means vibrating a system with 

frequency f and displacement amplitude b 

corresponds to a sinusoidal gravity modulation with 

the same frequency and acceleration amplitude ba> 2 

and vice versa (accordingly, hereafter the terms 

"gravity modulation", "periodic acceleration", 

"container vibration" and g-jitter will be used as 

synonyms). This also means that # a = borp T Md 3 rari j-, e 

va 

regarded as a variant of the classical Rayleigh number 
with the steady acceleration being replaced by the 
amplitude of the considered periodic acceleration. 

In general, the treatment of the problem is possible in 
terms of three independent nondimensional 
parameters only, where the first is the well-known 
Prandtl number (Pr) and the others are the 
nondimensional frequency (Q) and displacement (A), 
defined as: 



cod" 



A = b 



a 

J3 T AT 



(46) 



(47) 



where, obviously AQ 2 =PrRa a 



In the specific case of sufficiently small amplitudes 
(A«l) and sufficiently large frequencies (0»1) of the 
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vibrations, it is known (see the historical background 
reported in Lappa, 2010) that, for given Prandtl 
number, vibrational convection depends only on one 
relevant dimensionless parameter, the socalled 
vibrational Rayleigh number originally introduced by 
Gershuni (Gershuni and Zhukhovitskii, 1979): 



Ra v = 



(bco/3 T ATd) 2 (/3 T ATdY 



2va 



2va 



CO 



D. 2 A 2 



2Pr 



(48) 



where AT is the considered applied temperature 
gradient. 

Under the assumptions of small amplitudes (A«l) 
and large frequencies of the vibrations (Q»l), the 
Gershuni's formulation leads to a closed set of 
equations for the time-averaged quantities. The time- 
averaged continuity and energy equations remain 
formally unchanged (i.e. they correspond to eqs. (1), 
and (3), respectively, with t = t and v = v); the time- 
averaged momentum equation must be re-written as: 



a 



= -Vp-V 



vv 



+ Pr V 2 V + Pr Ra v [(vv • VT)n -w-Vvv] 

(49) 



where w appearing in the production term (the 
mathematical details related to the derivation of this 
term are given in Lappa, 2010) is an auxiliary potential 
function satisfying the equations: 



V-w = 

Vaw = VTa«->V 2 w = -Va 

y 



(vr An) 



(50a) 
(50b) 



//////////// 



/ / / / / 



/ / / / / ^/ / / / / / / 




//////////////////// 



FIG. 9 SKETCH OF LAYER OF INFINITE EXTENT SUBJECTED TO 
VIBRATIONS. 

Following Birikh (1990), the mathematical problem 
related to the derivation of an analytic solution for the 
conditions corresponding to Fig. 9 can be defined as 
follows. 

Assuming a generic plane-parallel flow solution in the 
form: 

u(y) 

(51) 



V 








w = 



(52) 
(53) 



T = x+f(y) 
w(y) 



the original partial differential equations for 

momentum, energy and the auxiliary potential 

function w given before can be reduced to a system of 

ordinary differential equations (here 9 is the angle 

between the shaking direction and the x-axis): 

d 3 u dw 
— + Ra v —cos(0) 
dy dy 

d 2 f 
dy 2 
dw 







- + u 







(54) 



(55) 



4f 



(56) 



= ^-cos(6>)-sin(6>) 
dy dy 

that combining eq. (54) with eq. (56) can be cast in 
compact form as: 



d 3 u 

d^£ 

dy 2 
with 



+ Ra v ^ cos 2 (6>) = Ra v cos(0) sin(0) 
dy 

+ u =0 



(57) 



(58) 



Ra v 



(ba/3 T yd 2 f 
2va 



where y is the rate of uniform 



temperature increase along the x axis and boundary 
conditions at y=±l/2: u=0 (no-slip) and d//dy=0 or f=0 
for adiabatic or conducting walls, respectively. 
Notably, the dimensionless characteristic number of 
the problem together with the components of the unit 
vector n appearing in these equations can be grouped 
in two parameters only: 



R 1 = Ra v cos(0)sm(0) 
R 2 =Ra v cos 2 (0) 



(59) 
(60) 



of which, the first one can be used to obtain conditions 
corresponding to the existence of the socalled "states 
of quasi-equilibrium" (no time-average flow) by 
simply setting it equal to zero (which leads to two 
possible cases of quasi equilibrium: the well-known 
condition of vibrations parallel to the temperature 
gradient (9=0°— >sin(9)=0-»Ri=0) and vibrations 
perpendicular to the temperature gradient 
(9=90°h>cos(9)=0^ Ri=0)). 

In general, the analytic form of the plane-parallel flow 
established in the layer for other values of 9 not 
corresponding to equilibrium (9^0°, 9^90°) must be 
determined by solution of the system (54)-(56) with the 
additional constraints: 



ju(y)dy 



and "r 2 . . , 

J w(y)dy = 

-l ; 2 



(61) 
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As illustrated in detail by Birikh (1990), this approach 
leads to the following mathematical expressions: 



u(y) 



16 



cosh(2£y) sin(2£y) sinh(!) cos(!)) 



! 3 (sinh(!)cosh(!) + sin(!)cos(!)) 

sinh(2gy) cos(2gy) cosh(!) sin(!) 
! 3 (sinh(!) cosh(!) + sin(!) cos(!))_ 



(62) 



2y- 



sinh(2£y) cos(2^y) sinh(!) cos(!) 
!(sinh(!) cosh(!) + sin(!) cos(!)) " 



cosh(2gy) sin(2gy) cosh(!) sin(!) 
!(sinh(!) cosh(!) + sin(!) cos(!)) 



(63) 



for the case with adiabatic walls 

cosh(2fo) sin(25Q sinh(!) cos(!) 
! 2 (sinh 2 (!) cos 2 (!) + cosh 2 (!) sin 2 (!)) 

sinh(2!y) cos(2!y) cosh(!) sin(!) 
! 2 (sinh 2 (!) cos 2 (!) + cosh 2 (!) sin 2 (!))_ 



«O0 = -— 

16 



(64) 



/OO = -|tan(0) 



2y- 



sinh(2!y) cos(2!y) sinh(!) cos(!)) 
(sinh 2 (!) cos 2 (!) + cosh 2 (!) sin 2 (!)) 

cosh(2!y)sin(2!y)cosh(!)sin(!) , 



(sinh 2 (!) cos 2 (!) + cosh 2 (!) sin 2 (!)) 

for the case with conducting walls 
where 

£ = ^[jto v cos 2 (60f /4 



(66) 



Modulated Buoyant Flows 

There have also been studies on the effect of time- 
modulated gravity (or vibrations) on systems with 
basic buoyant convection induced by vertical static 
(steady and uniform) gravity and horizontal 
temperature gradients 

The simplest model for this kind of flows is the 
canonical layer of infinite horizontal extent already 
examined before for pure buoyancy (the Hadley flow) 
and vibrational flow (pure thermo vibrational flow). 

By analogy with the analogous problem treated for the 
pure thermovibrational case, let us start the related 
discussion by observing that under the effect of high- 
frequency vibrations (in the framework of the 
averaged formulation discussed before) also this 
system admits solutions corresponding to quasi- 
equilibrium, i.e. states in which the mean velocity is 
zero. These states are made possible by a perfect 
balance of the static component of the body force and 
pressure gradient established in the liquid. The related 
mathematical (necessary) conditions of existence for 
the specific case considered here, i.e. vertical steady 



gravity and horizontal temperature gradients yr„ = i x , 
reduce to: 

Ra v cos(6») sin(0) - Ra = (67) 
When requisites for mechanical quasi-equilibrium are 
not satisfied time-average fluid motion arises. Such 
convective states are twodimensional and admit 
analytic solution in the form of plane-parallel flow. 
Following the general concepts introduced before for 
pure thermovibrational flows, these solutions can be 
represented as 



u(y) = - 



16 



cosh(2!y) sin(2!y) sinh(!) cos(!) 
! 3 (sinh(!) cosh(!) + sin(!) cos(!)) 



sinh(2!y) cos(2!y) cosh(!) sin(!) 
! 3 (sinh(!) cosh(!) + sin(!) cos(!)) 



(68) 



f(y)- 



2 R„ 



2y- 



sinh(2!y) cos(2!y) sinh(!) cos(!) 
!(sinh(!) cosh(!) + sin(!) cos(!)) " 



cosh(2!y) sin(2!y) cosh(!) sin(!) 
!(sinh(!)cosh(!) + sin(!)cos(!)) 

for the case with adiabatic walls 
R 



(69) 



u(y) : 



M 

16 



cosh(2!y) sin(2!y) sinh(!) cos(!)) 
! 2 (sinh 2 (!) cos 2 (!) + cosh 2 (!) sin 2 (!) ) 



sinh(2!y) cos(2!y) cosh(!) sin(!) 
! 2 (sinh 2 (!) cos 2 (!) + cosh 2 (!) sin 2 (!))_ 



(70) 



f(y)- 



2 R„ 



sinh(2!y) cos(2!y) sinh(!) cos(!) 
y (sinh 2 (!) cos 2 (!) + cosh 2 (!) sin 2 (!)) 



cosh(2!y) sin(2!y) cosh(!) sin(!) 
(sinh 2 (!) cos 2 (!) + cosh 2 (!) sin 2 (!))_ 



for the case with conducting walls (see Fig. 10) 
where 



(71) 



fi, = Ra v cos(<9)sin(60 - Ra 



R, 




(72) 

(73) 
(74) 



1 1 1 1 1 1 1 

-0.010 -0.005 0.000 0.005 
Temperature 



0.010 



a) 
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1 1 1 I I I 1 1 1 

-0.0002 0.0000 0.0002 

Temperature 5) 

FIG. 10 EXACT SOLUTION FOR THE CASE OF MIXED 
THERMOGRAVITATIONAL-THERMOVIBRATIONAL 
(AVERAGE) CONVECTION IN INFINITE LAYER WITH TOP 
AND BOTTOM SOLID CONDUCTING WALLS (9=0°, Ra=l): (a) 
VELOCITY PROFILE u(y); (b) TEMPERATURE PROFILE f(y). 

Mixed Marangoni/Thermovibrational Convection 

Like buoyancy flow, also the thermal Marangoni flow 
can be strongly affected by imposed vibrations. Such a 
case is considered in this section (surface-tension- 
induced flow interacting with vibrational effects in the 
absence of a static acceleration component, i.e. Ma^O, 
Ra m or Rav^O and Ra=0). 

Interesting results along these lines deserving 
attention are due to Suresh and Homsy (2001), who 
considered the infinite parallel Marangoni flow 
subjected to gravitational modulation at low 
frequencies (where the averaged model is not 
applicable) for the case in which finite-frequency 
vibrations are perpendicular to the layer (8=90°), (Fig. 
11; they assumed as base unmodulated Marangoni 
flow the popular "return flow solution" given by eq. 
(20) and employed a quasi-steady approach, in the 
limit of very low forcing frequency). 



Free surface 
u 




/////// / /// / // // // // // // // 



FIG. 11 SKETCH OF FLUID LAYER OF INFINITE EXTENT WITH 
UPPER FREE SURFACE SUBJECTED TO VIBRATIONS 
s(t) = bsin(cot)n WITH ARBITRARY DIRECTION. 



Following the general concepts introduced earlier, 
such solution (see Figs. 12) can be expressed as the 
superposition of two components, i.e. the steady 
Marangoni component proportional to Ma (crryd 2 /|ia) 
and a periodic vibrational component proportional to 
Ra m (ba) 2 Piyd 4 /va=AQ 2 /Pr, where y is the rate of 
uniform temperature increase along the x axis): 
u=Mag M (y) + Ra ia g B (y)exp(iQt) (75) 
T = x + Maf u (y) + RaJ B {y) exp(/Q0 (76) 

where, as already reported in the case of pure 
Marangoni flow, for adiabatic horizontal walls: 



g M (y) 



4 



3r + ,-i 



f M (y) = - 



j_ 

48 



3y + 2y ~2 y ~2 y -T6 



(77) 



(78) 



and, as analytically determined by Suresh and Homsy 
(2001): 



g B (y) = mc l exp(mr/) - mc 2 exp(-mr/) + c 3 - 
and: 

f B (y) = d x exp(«7/) + d 2 exp(-«7/) + d 3 exp(m?/) + 
d 4 exp(-mrj)+d 5 -rj/m 2 n 2 for Pr^l 
f B (y) = (d ] +d 3 r/)exp(mr/) + 
(d 2 +d 4 r/)exp(-mr/)+d s -r//m' 1 forPr=l 
with 



JL 

m 2 



71 = 



m = 



y + - 



c, = 



c, = 



Pr 

exp(-m)(m 2 /2 - 1) - m + 1 
2m 4 (-mcosh(m) + sinh(m)) 

exp(m)(-m /2 + 1) — m — 1 
2m 4 (-m cosh(m) + sinh(m)) 
m(c 2 - q) 



c 4 =-(c 1 +c 2 ) 



(79) 

(80a) 
(80b) 

(81a) 

(81b) 

(81c) 
(82a) 

(82b) 

(82c) 
(82d) 



if Pr^l 



1 



2«sinh(«) In —m 



m 



r[(c, + c 2 ) exp(-n) - c, exp(m) ■ 



, si 1 - exp(-n) 
c 2 exp(-m)J + fj-^- 



(83a) 



m 



2nsinh(n) [n 2 - m 

, . n 1 - exp(n) 
■c 2 exp(-m)J+ 



[(q + c 2 ) exp(«) - Cj exp(m) ■ 



(83b) 
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d 3 = - 



d, = 



5 n 2 



if Pr=l 



(83c) 
(83d) 
(83e) 



1 



[2 sinh(OT)c, + m{c l exp(m) ■ 



4msinh(m) 

2(1 - exp(-m)) 



c 2 exp(-m)) + 
1 



Am sinh(m) 

2(1 - exp(m)) 



[- 2sinh(m)c 2 + m(c, exp(m) - 



c 2 exp(-m)) + 



(84a) 

(84b) 
(84c) 
(84d) 
(84e) 



with the horizontal boundaries located at y=±l/2. 




0.50- 
0.40- 
0.30- 

0.20 
0.10-I 

>» 0.00 r 
•0.10- 
-0.20 
-0.30-3 
-0.40 
-0.50 



-0 40 -t) 20 0.00 

Velocity 

Ii. I I . 



0.20 




a) 



I I I I I I I J i I I I I I I I I I I I I I I I I 

-0.01 000 0.01 0.02 0.03 

Temperature b) 

FIGS. 12 EXACT SOLUTION FOR THE CASE OF MIXED 

MARANGONI-THERMO VIBRATIONAL CONVECTION IN 

INFINITE LAYER WITH ADIABATIC BOUNDARIES (Pr=l, Ma=l, 

6=90°, Ra„=5, Q=l). 



Solutions in the Presence of Static and 
Uniform Magnetic Fields 

In many material processing techniques, a magnetic 
field is commonly used to control the liquid flows. Its 
action can lead to the braking of the flow (i.e. the 
reduction of the rate of convective transport) or to the 
damping of possible oscillatory convective instabilities. 
For these reasons, a number of theoretical studies have 
appeared during recent years concerning the effect of 
magnetic fields on the basic flow motion in several 
geometrical models of semiconductor growth 
techniques. 

The present section illustrates the effect of a constant 
magnetic field (static and uniform) on the 
fundamental solutions for gravitational and 
Marangoni convection introduced in earlier sections. 

Physical Principles and Governing Equations 

The motion of an electrically conducting melt under a 
magnetic field induces electric currents. Lorentz forces, 
resulting from the interaction between the electric 
currents and the magnetic field, affect the flow. 

In the following, as stated in the title of the present 
subsection, the "physics" of the problem and the 
introduction of the corresponding model equations are 
considered for the simple and representative case of 
constant magnetic fields. 

A uniform magnetic field with magnetic flux density 
(also referred to as magnetic induction) Bo generates a 
damping Lorentz force through the aforementioned 
electric currents induced by the motion across the 
magnetic field that in dimensional form can be 
expressed as F L =j f AB where j is the electrical 
current density. Like other body forces (e.g., the 
buoyancy forces), an account for this force can be 
yielded simply adding a relevant term to the 
momentum equation. 

Following this approach, scaling the magnetic flux 
density with Bo and the electric current density j 

with rjeVrefBo, the momentum equation including the 
Lorentz force, can be written in nondimensional form 
(assuming as usual Vref=a/d) and in the absence of 
phase transitions as: 



a 



-Yl> -V-[yy]+PrV 2 V-Pr RaTi g +PrHa 2 {j_ f a/ b J 



(85) 



where Ha is the socalled Hartmann number 
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Ha = B o d\^-j 



(86) 



cje is the electrical conductivity and iBo the unit vector 
in the direction of Bo. 

The nondimensional electric current density / can be 
expressed via the Ohm's law for a moving fluid as: 

l f =E + VAi Bo (87) 

where E is the electric field normalized by VrefBo. Since 
(see Lappa, 2010), in many cases of technological 
interest) the unsteady induced field is negligible, in 
particular, the electric field can be written as the 
gradient of an electric potential: 

£ = -Vd) c (88) 
The conservation of the electric current density gives: 

Y-l f =0 (89) 
which combined with eq. (88) leads to a Poisson 
equation for the electric potential: 

V 2 O e =V.(VA/J = / Bo .y>V (90) 
Finally, the momentum equation can be rewritten as: 

dV 



= -Vp - V • \VV] + Pr V 2 V - Pr RaTi 



(91) 



a 

+ Pr//a 2 (-YO e +yA/ B jA/ Bo 

Buoyant Convection with a Magnetic Field 

The simplest model in such a context is the flow in an 
infinite planar liquid metal layer confined between 
two horizontal solid walls driven by a horizontal 
temperature gradient (the canonical infinite Hadley 
flow whose general properties have been treated in a 
previous section in the absence of magnetic fields). 

Notably, for electrically insulating horizontal walls 
such a flow admits analytical expression even in the 
cases in which it is subjected to a magnetic field with 
direction vertical or parallel to the applied 
temperature gradient, or even directed spanwise to the 
basic flow (see Fig. 13). 




FIG. 13 SKETCH OF LIQUID LAYER UNDER THE EFFECT OF A 
HORIZONTAL TEMPERATURE GRADIENT AND A MAGNETIC 



FIELD WITH VERTICAL, STREAMWISE OR SPANWISE 
DIRECTION. 

In the case of vertical field, the analytic solution (see, 
e.g., Kaddeche et al., 2003) reads: 



u = - 



Ra 



Ha 1 



sinh(Hay) 
2sinh(/fo/2) 



-y 



x + - 



Ra 

Ha 2 



1 sinhCggy) y | 
2Ha 2 sinh(/fa/2) 6 



(92) 



(93) 



with the horizontal boundaries located at y=±l/2 and 
the parameter r\ given by: 

1 cosh(//a/2) 
8 2Hasmh(Ha/2) 



1 



(94a) 



for adiabatic boundaries and 



1 



1 



24 Ha 2 



(94b) 



for conducting boundaries. 



If a horizontal magnetic field parallel to the applied 
temperature gradient (magnetic field applied along 
the x-direction) is considered, there is no direct effect 
of the field on the parallel flow in the layer as the 
velocity is parallel to the field direction (vAi Bo =o), 
which (taking into account eqs. (87)-(90)) leads to the 
vanishing of the production term in eq. (91). The basic 
flow is, hence, the flow without magnetic field given 
byeq.(18). 

For a magnetic field still horizontal but directed 
spanwise to the basic flow (i.e. a magnetic field 
applied along the z-direction), the basic parallel flow 
in the layer is still the flow without magnetic field. 

A justification for such a behavior is not 
straightforward as one would assume, and deserves 
some specific additional insights. 

First of all, it has to be noted that as long as we 
consider the basic flow as being horizontally 
homogeneous, the electric field induced by the liquid 
motion is uniform and directed perpendicular to the 
plane of the layer. Zero circulation of this field (related 
to the fundamental irrotational nature of electric fields 
when steady magnetic fields are considered) implies 
that no electric currents closing within the liquid layer 
can be induced, but the electric impermeability of the 
horizontal boundaries also precludes the possibility of 
any electric current passing normally through the 
layer; thereby, insulating horizontal walls lead to a 
separation of electric charges over the depth of the 
layer, so giving rise to an electrostatic field which 
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cancels that induced by the liquid motion. As a result, 
there is no electric current induced by a horizontally 
uniform basic flow in a coplanar magnetic field and, 
therefore, there is no direct influence of the magnetic 
field on the basic flow. 

Marangoni Flow with a Magnetic Field 

Notably, like the Hadley flow discussed before, also 
the infinite surface-tension-driven flow admits 
analytical expression under the effect of a magnetic 
field. 

In particular, if the magnetic field is parallel to the free 
surface, i.e. satisfies Biy=0 (hereafter referred to as 
coplanar field), then, following the same arguments 
already provided for the infinite Hadley flow, it can be 
concluded that the field has no influence on the basic 
flow, and consequently the parallel flow remains the 
same as without it (this is exactly the return flow given 
byeq.(20)). 

By contrast, if the magnetic field is perpendicular to 
the free surface the analytic solution (neglecting 
exponentially small terms for strong magnetic field 
Ha»l) becomes (Priede et al., 1995): 



-Mai— — exp 
[Ha 



y Ha 



Ha 



-exp 



-\\+y\Ila 



Ha' 



(95) 



T = x-Ma{- 



y ( y + i 



Ha 1 



Ha 



-exp 



-\)--y\Ha 



y i 

-^T + rexp 

Ha 3 Ha 4 



+ y \Ha 



(96) 



where, as usual, the horizontal boundaries are 
assumed to be located at y=±l/2 and Ma= crryd 2 /|ia (y 
being the rate of uniform temperature increase along 
the x axis and d the depth of the layer), Ha=Bod 

(Ge/pv) 1 ' 2 . 
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